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We describe hard exclusive processes involving a transversally polarized 
p meson in the twist 3 approximation, in a framework based on the Taylor 
expansion of the amplitude around the dominant light-cone direction. 
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1. Introduction 

The last 15 years have witnessed a tremendous development of the QCD 
understanding of hard exclusive processes. The basis of this progress has 
been the derivation of factorization proofs [1] for exclusive amplitudes in 
various generalized Bjorken kinematical regimes, at the leading twist level. 
Schematically, the amplitude of a process governed by a large energy scale 
Q is written as the convolution of a perturbatively calculable subprocess 
amplitude and a few hadronic matrix elements of light cone operators. For 
instance, near forward hard electroproduction of a p meson depends on the 
H(x, £, t) and E(x, £, t) GPDs and on the p distribution amplitudes (DA) 
<S> p (z). 
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To extend these factorization procedures beyond the leading twist is 
mandatory in a number of cases including the phenomenologically impor- 
tant instance of the production of transversally polarized vector mesons. 
The understanding of the quark-gluon structure of a vector meson is an 
important task of hadronic physics if one cares about studying confinement 
dynamics. This quark gluon structure may be described by distribution 
amplitudes which have been discussed in great detail [2]. 

Contrarily to the longitudinally polarized case, the leading twist 2 DA of 
the transversally polarized vector mesons is chirally-odd, a property which 
opens the way to its use for the study of chirally-odd GPDs [3j (although 
it is often decoupled because of the vanishing of hard amplitudes [I] ). One 
thus needs to calculate amplitudes at the twist 3 level [5j [6] and this is the 
source of quite a number of subtle problems, as already discussed in [I]. 
Two issues are worth emphasizing. Firstly, new non-perturbative objects 
enter the discussion, as for instance matrix elements of operators containing 
a transverse derivative. Secondly, the end-point singular behaviour requires 
careful treatment of kinematical region with soft partons. 

Both difficulties should be carefully addressed in turn. Here, we shall 
not discuss the second one, which may hopefully be treated from a modified 
collinear factorization point of view [5], including partonic transverse de- 
grees of freedom and subsequent Sudakov resummation techniques. We plan 
to return to this question later. The first difficulty has been overviewed by 
most authors, leading sometimes to the use of gauge dependent expressions 
for the amplitudes. Indeed preserving gauge invariance in a factorization 
procedure is not guaranteed, as soon as one allows transverse motion of the 
partons at some intermediate step of the calculation. To preserve the QCD 
gauge invariance, a subtle interplay of various contributing amplitudes must 
conspire. For instance, in the case we will detail below, QCD gauge invari- 
ance is translated to the statement that the impact factor should vanish 
when the virtuality of an exchanged t— channel gluon vanishes. We will 
show that this is indeed the case provided two and three particle correlators 
are simultaneously taken into account. 

Although some courageous attempts have tried to develop the phe- 
nomenology of exclusive vector meson production [SJ [H] a careful analysis of 
the specificities of the transversally polarized case is still lacking. The ex- 
isting data from HERA (both at high energy with HI and ZEUS detectors, 
and at lower energy with HERMES) on the one side, the forthcoming data 
from JLab and Compass on the other side, deserve a complete study which 
we aim to perform within the twist 3 approximation. As a first step [9], we 
show how to use properly the two and three particle distribution amplitudes 
within a factorized approach. 
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2. Calculation of the 7* — ► pT impact factor 

As the simplest example of a hard exclusive reaction involving a transver- 
ally polarized vector meson, let us consider the high energy process 

7*(?) + 7V) -» Pr(pi) + p(P2) (1) 

where the photons are highly virtual : q 2 ,q' 2 = —Q 2 , —Q' 2 3> ^-qcd > an d 
the Mandelstam variables obey the conditions s S> Q 2 , Q' 2 , — i- Neglecting 
meson masses, one considers for reaction ((!]) the vector meson momenta as 
the light cone vectors (2pi • p 2 = s); the "plus" light cone direction being 
directed along p\ and the "minus" light cone direction along p 2 with vector n 
defined as P2/(pi ~P2)- In this Sudakov basis, transverse euclidian momenta 
are denoted with underlined letters. The virtual photon momentum q reads 
q = pi — n. The impact representation of the scattering amplitude A4 
for the reaction ([!]) is 



.-> 2 / 5+ioo 



(2^/ -f^r_-k)\ ^f(-g, -r_ + k>)j — (-) Cfc&'.E) 

5—ioo 



where G w is the 4-gluons Green function which obeys the BFKL equation. 
G w reduces to l/u)5{k — !±)k 2 /(r ~ t} 2 within Born approximation. We 
focus here on the 7* — ► p impact factor $ which is the integral of the k- 
channel discontinuity of the S-matrix element <S^ T of the subprocess 

g(h,ei) + 7*(g) -> g(k 2 ,e 2 ) + Pt(pi) 

d>7*-P(^ £ _ k ) = e 7V _ / _ DigCK 5 7* L _ k): ( 2 ) 

where k = (g + A;i) 2 . Considering the forward limit for simplicity, the gluon 
momenta read 

7 /£ -|- Q 2 ~\~ k_ K -\- k 9 ,0 T 9 ?2 

«i = P2 + »T) k 2 = P2 + kx, k x = k 2 = k T = -k . 

s s 

The impact factor $ can be calculated within the collinear factorization. 
It is a convolution of perturbatively calculable hard-scattering amplitudes 
H and soft correlators S involving relevant p-meson distribution amplitudes, 
symbolically written as 

$ = j dH --- tr[H(l---)S(l---)} , (3) 

and which involves loop integrations formed by the n partons (n > 2) enter- 
ing soft correlators S. Working within the light-cone collinear factorization 
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Fig. 1. 2- and 3-body correlators attached to a hard scattering amplitude in the 
example of the 7* — ► p impact factor. 



framework, let us first derive the contribution of the diagrams with the 
quark-antiquark correlators. The basic tool is to decompose any hard coef- 
ficient function H(l) around a dominant "plus" direction: 



dl 



a 



(e-yp) a + ..., (4) 

C=yp 



with (l — yp) a ~ (-a U P to twist 3 accuracy. Here and subsequently a 
denotes transverse components. This collinear expansion of the hard part 
should be compared to the expansion of matrix elements in powers of space- 
time separation (x 2 ) n in the coordinate approach. 

The I 1 - dependence of the hard part looks first as an excursion out of 
the collinear framework. The factorized expression is obtained as the result 
of an integration by parts which replaces £^ °y acting on the soft corre- 
lator. This leads to new operators O 1 - which contain transverse derivatives, 
such as 'tpd^tp, and thus to the necessity of considering additional correla- 
tors ^ ± (l). This procedure leads to the factorization of the amplitude in 
momentum space. Factorization in the Dirac space can be achieved by a 
Fierz decomposition. Thus, the amplitude takes the factorized form, see 
Fig. 1: 



1 1 

A = J dytT[H(y)T] <D r (y) + J dy 1 dy 2 tv[H^y 1 ,y 2 )T] $J(yi,jft). (5) 



The parameterizations of the vacuum-to-/j-meson matrix elements needed 
in the above factorization procedure up to the twist-3 order in the axial 
(light-like) gauge n ■ A = 0, which we will use in order to get rid of Wilson 
line, can be written as (here z = An) 

(p(p)M*)7^(0)|0> = m p / p [^ 1 (y)(e*-n)p /1 + ^ 3 (y)e; T ], 
<P(P) 1^)757^(0)10) = J m p f p 

and 

(p(p)Mz)W 5^(0)|0> ^ m p f p ^(y)p,e* a T , 
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1^0)757/^ ip(0)\0) = rn p f p iip^(y)p fl e a xi38e*x T p l 3ns, (7) 

where = denotes the Fourier transformation dye^ typ ' z ^ . The matrix ele- 
ments of quark-gluon nonlocal operators can be parameterized as 

(pipmz^gAKz^miO) = m p fV p B( yi ,y 2 )p,ef, (8) 

(p(jP)|$(2l)757MS^ate)^(0)|0) = m p fz p i D(y 1 ,y 2 ) p^axps ef Pp n s , 

11 T 
where / dy x J dy 2 e^f^-K^-wJf^] i s denoted by = . The light-cone frac- 
o o 

tions of the quark, anti-quark and gluon are y\, 1 — y 2 and y 2 — yi ■ 

The correlators introduced above are not independent, thanks to the 
QCD equations of motion for the field operators entering them (see, for 
example, [6]). In the simplest case of fermionic fields, they follow from the 
vanishing matrix elements {iD(0)ip(0) ip{z)) = and {ip(0) iD(z)vp(z)) = 
due to the Dirac equation. Denoting C^'^ = f^ p A / f p , we obtain 

m fsivi) + m PA(yi) + vliyi) + <^a(j/i) = 
i 



J dy 2 [C 3 y B( V1 , y 2 ) + C 3 A D( yi , y 2 )] (9) 



and 



yi <pz(yx) - y\ VA{yi) - (yi) + pa(vi] 
i 



J dy 2 B(y 2 , y x ) + C 3 A D(y 2 , Vl )] . (10) 



A crucial point of this approach is that scattering amplitudes do not 
depend on the specific choice of the vector n which fixes the light-cone 
direction. As observed for inclusive structure functions [10] and for DVCS 
|llj , this n independence condition leads to non-trivial relations at the twist 
three level. In the case of exclusive processes [6], this constrains the non- 
perturbative correlators entering the factorized amplitude. The polarization 
vector for transverse p which enters in the parametrization of twist 3 corre- 
lators depends on n according to 

=e /1 -p fl e -n (11) 

i.e. it is defined with respect to the light-cone vector n. The n— independence 
condition can thus be written as 

dA , d d . d 

— = 0, where — = — - + e*— -. 12 

dn» dn» dn^ M d(e* ■ n) K ' 



G 
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Fig. 2. Reduction of 3-body correlators to 2-body correlators through Ward iden- 
tity. 

The appearance of the total derivative may be interpreted as a (vector) 
analog to the renormalization group (RG) invariance when the dependences 
on the renormalization parameter coming from various sources cancel. One 
can view this as a RG-like flow in the space of light-cone directions of con- 
tributions to the amplitude where the polarization vector plays the role of 
a beta function. This condition expressed at the level of the full amplitude 
can be reduced to a set of conditions involving only the soft correlators. 
The general strategy relies on the use of the collinear Ward identities which 
relate firstly amplitudes with different number of legs and secondly higher 
order coefficients in the Taylor expansion ^) to lower order ones. The most 
involved use of these identities occurs in the case of the 3-body correlator. 
In the case of vector correlator (|5J), due to (fTTj) the dependency on n enters 
linearly and only through the scalar product e* • n . Thus, the action on the 
amplitude of the derivative d/dn defined in (fT2j) can be extracted by the 
replacement e* — > p a , and after using the Ward identity, it reads 

{y 1 -y 2 )tv[H p {y 1 ,y 2 )pPl>]=tT[H(yi)f [ -tY{H(y 2 )i>) , 

as illustrated by Figj2j One can apply similar tricks to the 2-body correla- 
tors with transverse derivative whose contribution can be viewed as 3-body 
processes with vanishing gluon momentum. This finally leads to the equa- 
tions 

d 1 
-r-Vi{yx) = -(Piti/i) +^3(yi) - (3 [ dm (B(yi,y 2 ) + B(y 2 ,yi)) , 
dyi J 2/2-2/1 

1 

^ T A (yi) = Mvi) - (£ I -^-(D( yi ,y 2 ) + D( y2 , yi )) . (13) 

dyi J y 2 - yi 



We emphasize that the equations obtained above for the special case of 
impact factor are universal i.e. do not refer to any specific hard process. 
The equations of motion together with those expressing the n-independence 
are four equations which constrain seven DAs (</?i, ip3, <pi, ifA, <p^, B, D). 
Thus, any hard process involving the exclusive production of a p can be 
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expressed, at twist 3, in terms of the three independent DAs which we 
choose as (p±, B and D. The solution of Eqs. (f9l [T0l[T3j) can be generically 
written as the sum of the solution in the Wandzura-Wilczek approximation, 
i.e. with vanishing B and D DAs, and the solution describing the genuine 
twist 3 part. The Wandzura-Wilczek parts take the forms 



,„TWWf„,\ 
"Pi/ A W 



dv 
v 

y 



. . f dv . , 
ipi(v) ± / — tpi(v) 



-!J I ^r<Pi(v) ± y 



dv 



-tpx[v) 



(14) 
(15) 



The remaining genuine twist 3 DAs read: 



J y 2 -u J y 2 ~ U 



dy 2 



1 f du 



/f Jdy 2 ±(CYB + tiD)(u,y 2 ) 



(16) 



/ -^-(tf* + C?D)(u, y 2 ) - [ -^{$B + CiD)(y 2 , 
J V2-u J y 2 - u 

B{yi, y 2 ) 



u 



y y i 

¥i gen (y) = / duipf r \u) -(3 dyi / dy 2 











v 



V2 - y\ 



(17) 



while the corresponding expressions for ip A en (y) and f\ 9en {y) are obtained 
by the substitutions: 



qeni \ S3 D S3 u 



<p$T(.y) , 



Tqerii \ S3 " " S3 ^ l qen/ \ 

<pa (y) < — > (y) 



(18) 
(19) 



In summary, we describe consistently exclusive processes at higher twist, 
in a way which explicitly preserves gauge invariance. The n-independency 
condition generalized up to the dynamical twist 3 plays a crucial role for the 
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consistency of this approach with the studies of DAs performed in [2]. The 
present framework opens the way to a systematic and consistent treatment 
of hard exclusive processes. This does not preclude the solution of the end- 
point singularity problem 0(8] which requires a dedicated treatment. 
We thank O. V. Teryaev for discussions. This work is supported by the 
Polish Grant N202 249235, the French-Polish scientific agreement Polonium, 
the grant ANR-06-JCJC-0084, the RFBR (grant 09-02-01149 and 07-02- 
91557, 08-02-00334, 08-02-00896). 
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